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Abstract 

Nonlinear convex relaxations of the power flow equations and, in particular, the Semi-Definite Pro¬ 
gramming (SDP), Convex Quadratic (QC), and Second-Order Cone (SOC) relaxations, have attracted 
significant interest in recent years. Thus far, little attention has been given to simpler linear relaxations 
of the power flow equations, which may bring significant performance gains at the cost of model accu¬ 
racy. To fill the gap, this paper develops two intuitive linear relaxations of the power flow equations, one 
based on classic network flow models (NF) and another inspired by copper plate approximations (CP). 
Theoretical results show that the proposed NF model is a relaxation of the established nonlinear SOC 
model and the CP model is a relaxation of the NF model. Consequently, considering the linear NF and 
CP relaxations alongside the established nonlinear relaxations (SDP, QC, SOC) provides a rich variety 
of tradeoffs between the relaxation accuracy and performance. 


Nomenclature 


N - The set of nodes in the network 
E - The set of from edges in the network 
E r - The set of to edges in the network 
i - imaginary number constant 
I - AC current 
S = p + iq - AC power 

V — v/L9 - AC voltage 

Z = r + ix - Line impedance 

Y = g + ib - Line admittance 

T = tZ6 t - Transformer properties 
Y a — g a + ib s - Bus shunt admittance 
W - Product of two AC voltages 
L - Current magnitude squared, |/| 2 
b c - Line charging 


s u - Line apparent power thermal limit 
9 a - Phase angle difference limit 
S d = p d + iq d - AC power demand 
S 9 — p 3 + iq 9 - AC power generation 
Co, ci, C 2 - Generation cost coefficients 
K(-) - Real component of a complex number 
Sy(-) - Imaginary component of a complex number 
(•)* - Conjugate of a complex number 
■ | - Magnitude of a complex number, Z 2 -norm 
Z - Angle of a complex number 
x u - Upper bound of x 
x l - Lower bound of x 
x - Convex envelope of x 
x - A constant value 
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1 Introduction 


Nonlinear convex relaxations of the power flow equations have attracted significant interest in recent 
years. They include the Semi-Definite Programming (SDP) [2], Second-Order Cone (SOC) [10], Convex- 
DistFlow (CDF) [6], and the recent Quadratic Convex (QC) [8] and Moment-Based [16, 17] relaxations. 
Much of the excitement underlying this line of research comes from the fact that the SDP relaxation 
has shown to be tight on a variety of case studies [12], opening a new avenue for accurate, reliable, and 
efficient solutions to a variety of power system applications. Indeed, industrial-strength optimization tools 
(e.g., Gurobi [7], cplex [9], Mosek [23]) are now available to solve various classes of convex optimization 
problems. 

Thus far, study of power flow relaxations has focused primarily on these nonlinear methods in the 
interest of the highest possible model accuracy at any cost. Subsequently, little attention has been given to 
simpler linear relaxations [22, 4], which are less accurate but have significant performance and scalability 
benefits. To fill that gap, this paper develops two intuitive linear relaxations of the power flow equations, 
one based on classic network flow models [1] and another inspired by copper plate approximations. The 
main contributions of this work can be summarized as: 

1. Developing two intuitive linear relaxations for AC transmission system optimization, the network 
flow (NF) and copper plate (CP) relaxations. 

2. Proving that NF is a linear relaxation of the established nonlinear SOC relaxation and that CP is 
a relaxation of NF. 

It is important to emphasize that, contrary to the models presented herein, the established DC power 
flow model [21] is a linear approximation of the AC power flow, not a relaxation. As a consequence, it 
cannot be used for providing quality guarantees (i.e. lower bounds) on the underlying AC power flow 
model. 

The rest of the paper is organized as follows. Section 2 reviews the extended formulation of non- 
convex AC power flow feasibility problem for transmission systems (AC-E-PF). Section 3 develops the 
relaxations of (AC-E-PF), it begins with a review of the SOC relaxation and then proposes the NF and 
CP relaxations. Section 4 reviews the related works on linear relaxations and illustrates how all of the 
relaxations considered perform on a well known 3-bus example. Section 5 concludes the paper. 


2 Extended AC Power Flow for Transmission Systems 

In the interest of clarity, AC Power Flows, and their relaxations, are most often presented on the simplest 
version of the AC power flow equations. However, transmission system test cases include additional 
parameters such as bus shunts (Y B ), fine charging (fo c ), and transformers (T), which complicate the AC 
power flow equations significantly. This section presents an extended AC power flow problem including 
all of these additional parameters. 

A power network is composed of a variety of components such as buses, fines, generators, and loads. 
The network can be interpreted as a graph (N, E ) where the set of buses N represent the nodes and the 
set of lines E represent the edges. Note that E is an undirected set of edges, however each edge (i, j) e E 
is assigned a from side (i,j) and a to side (j,i), arbitrarily. These two sides are critically important as 
power is lost as it flows from one side to another. Lastly, to break numerical symmetries in the model 
and to allow easy comparison of solutions, a reference node r e N is also specified. 

The AC power flow equations are based on complex quantities for current I, voltage V, admittance 
Y, transformers T, and power S, which are linked by the physical properties of Kirchhoff’s Current Law 
(KCL), 


Ohm’s Law, i.e., 


7? — if — YiVi = Y, Aj + £ 7ji Vi e N 

(i,j)EE ( j,i)EE 

Iij = (V, • i'Y) ^ y,,V, V(*,j), (. hi) ■- E 


(1) 


( 2 ) 
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and the definition of AC power, i.e., 


Sa = V(i, j) e E 

± ij 

Combining these three properties yields the extended AC Power Flow equations, 
Sf-S?-(1T)*|^| 2 = £ Sy+ £ SjiVieN 

(i,j)EE U,i)EE 

b f \ IWI2 T/.W* 


&,■ = V* - i- 


l^l 2 * T)* 


% = ^ \Vj\ 2 - Y, 


:V*Y 


(i,j) e E 


(3) 


(4a) 

(4b) 

(4c) 


Observe that ^ over (i,j) e E collects the edges oriented in the from direction and ^ over (j,i) e E 
collects the edges oriented in the to direction around bus i e N. These non-convex nonlinear equations 
define how power flows in the network and are a core building block in many power system applications. 
Additionally, many applications require various operational side constraints on the flow of power. We 
now review some of the most significant ones. 


Generator Capabilities AC generators have limitations on the amount of active and reactive power 
they can produce S 9 , which is characterized by a generation capability curve [11]. Such curves typically 
define nonlinear convex regions which are often approximated by boxes in AC transmission system test 
cases, i.e., 

sf SS S? < Sf u Vie N (5) 

Line Thermal Limit AC power lines have thermal limits [11] to prevent lines from sagging and 
automatic protection devices from activating. These limits are typically given in Volt Amp units and 
constrain the apparent power flows on the lines, i.e., 

V(i,i),(i,*)e£ (6) 

Bus Voltage Limits Voltages in AC power systems should not vary too far (typically +10%) from 
some nominal base value [11]. This is accomplished by putting bounds on the voltage magnitudes, i.e., 

v\ < [Vi[ Vie N (7) 

Phase Angle Differences Small phase angle differences are also a design imperative in AC power 
systems [11] and it has been suggested that phase angle differences are typically less than 10 degrees in 
practice [19]. These constraints have not typically been incorporated in AC transmission test cases [25]. 
However, recent work [5, 8] have observed that incorporating Phase Angle Difference (PAD) constraints, 
i.e., 

- 9* < ^ 0* V(i,j ) e E (8) 

is useful in the convexification of the AC power flow equations. For simplicity, this paper assumes that 
the phase angle difference bounds are symmetrical and within the range (— tt/2, 7t/ 2), i.e., 

\ (i,j)eE 

but the results presented here can be extended to more general cases. Observe also that the PAD 
constraints (8) can be implemented as a linear relation of the real and imaginary components of V, V* 
[15], i.e., 

tan(-0£)5R (ViV*) ^ 3 (ViV*) ^ tan(6»^)K (ViV*) V(i,j) e E (9) 
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Model 1 The Extended AC Power Flow Feasibility Problem (AC-E-PF) 


variables: 

Sfe(Sf l ,Sf u ) Vi e N 

Vie(y?,vn VieN 

S.jC iS'^.S") V{i,j),(j,i)eE 

subject to: 

ZF r = 0 

v\ ^ \V t \ VieN 

S?-St-(Yfi* |V*| 2 = ^ S ij+ ^ Sji Vie N 

(i,j)eE (j,i)zE 


II 

to 

\Y* 

1*7 

o* 

|£: n 

, N 2 

ViV; 

1 3 

2 ) 

' l^l 2 

*7 rr\* 

ij 

ii 

(y* 

1 *7 

S>. 

1 O'* 

to 

1 

y/•* v l v 3 
**7 r r. . 

J -13 


ISiiKs" V(i,j),(j,i)eE 

<>* < /(V.Vf) * 0* V(i,j)eE 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 

(13f) 

(13g) 


and that equation (4b) can be used to express these in terms of the S variables as follows, 

bf 


ViVj = Z,- T* 


Y' — i -44 

*3 2 


\Vi \ 2 

\Ti , I 2 * y 


(i,j) e E 


( 10 ) 


These equations combined with (9) implement the PAD constraints in terms of the V and S variables as 
follows, 


tan(-fl£)H (z*T* (jV t * 
tan(0£)H (zfjr* ( ( F* 










V(i,j)eE 

( 11 ) 

V(i,j)eE 

( 12 ) 


The usefulness of these alternate formulations will be apparent later in the paper. 

Other Constraints Other line flow constraints have been proposed, such as, active power limits 
and voltage difference limits [12, 15]. However, we do not consider them here since, to the best of our 
knowledge, test cases incorporating these constraints are not readily available. 

The Extended AC Power Flow Feasibility Problem Combining all of these constraints yields 
the Extended AC Power Flow Feasibility presented in Model 1 (AC-E-PF). The operational constraints 
bus voltage and generator output are captured by the variable bounds. Constraint (13a) sets the reference 
angle, to eliminate numerical symmetries. Constraint (13b) capture the bus voltage limits. Constraints 
(13c) capture KCL and constraints (13d)-(13e) capture Ohm’s Law. Finally constraints (13f) and (13g) 
enforce the line flow and phase angle difference limits respectively. Notice that this is a non-convex 
nonlinear satisfaction problem due to the product of voltage variables, V t V* and is NP-Hard [24, 13]. 
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Noting that bounds on the voltage variables and edge power flow variables (i.e. Vi e (V , 1 , ),S i3 e 

(iS-j. )) are not always specified, we observe that constraints (13b), (13f) respectively imply reasonable 


bounds on each of these variables. 

Lemma 2.1. 5“ = , Slj = — (s“- + is“ ) are valid bounds in (AC-E-PF). 

Proof. First, observe that (13f) is equivalent to, 

K ( 5 tf ) 2 + 9 f ( Sij ) 2 <(««) 2 ( 14a ) 

solving for K(SV,) results in, 

K(^) 2 <(4) 2 -3(^) 2 (15a) 

- VK ) 2 -®(^) 2 < K < V ( fl y ) 2 - 3 ( s y ) 2 (15b) 

Noticing that (s“) 2 ,Sy(-) 2 are both positive, we can deduce the following bounds, 

— < 5R (Sij) A sTj (16a) 

A similar argument holds for the bounds of Sy (Sij), demonstrating the result. □ 

Lemma 2.2. V“ = v? + ivf, V'' :j = -(u“ + ivf) are valid bounds in (AC-E-PF). 

Proof. First, observe that the upper bound constraint in (13b) is equivalent to, 

K(T) 2 +S(T) 2 (t>“) 2 (17a) 

The result follows similarly to the previous one. □ 


Also observe that the lower bound constraint in (13b) is non-convex and cannot, in general, be incorpo¬ 
rated into the bounds on V. 


3 Three Relaxations of the Extended AC Power Flow 

This section develops three successive relaxations of (AC-E-PF). It begins with the established Second- 
Order Cone (SOC) relaxation [10]. The SOC relaxation is then further relaxed in to a Network Flow 
model (NF), similar to those traditionally studied in operations research. Lastly, the NF model is then 
further relaxed into a Copper Plate model (CP), which ignores all of the network aspects and simply 
states that the supply must be at least as large as the demand. 


3.1 The Second-Order Cone Relaxation (SOC) 

The SOC relaxation was first proposed in [10] and utilizes two key insights. First, by lifting the product 
of voltage variables V, V* into a higher dimensional space (i.e. the FF-space), 

Wi = \Vi \ 2 isN (18a) 

Wij = ViV* V(i,j)eE (18b) 


a linear relaxation of (AC-E-PF) is obtained. Second, a relaxation of the absolute square of the voltage 
products is developed to strengthen this FF-space relaxation, as follows, 


T-V'A ( 

T* 3 \ T* 3 

*3 / \ *3 


rn* V 3 I 1 rp . V J 


(19a) 


l ^*| 2 _ N 2 | V .|2 

\ T *\ 2 IT*I 2 21 

\ViV *\ 2 = \Vifm 2 
\Wijl 2 = WAV, V(i,j) e E 
\Wij \ 2 V WiWj V(i, j) e E 


(19b) 

(19c) 

(19d) 

(19e) 
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Model 2 The SOC Relaxation of Extended AC Power Flow (SOC-E-PF) 


variables: 


S?e(Sf l ,Sf u ) Vi e N 
Wi e ((vl) 2 , «) 2 ) Vi e N 
W tj V(i,j) g E 


subject to: 


v(», 

i ) e 

E 


■ = E 

Sij + 


(>j)eE 

(j- 

Y *-i b A) 

Wi 

V* Wij 

. y 2 ) 

\Tij\ 2 

ij rr\* 

ij 

b c ), \ 


w* 

f) 

Wj - Y* 

(i 

T ■ v 

J-ij 




\Wijl 2 < WiWj V(i,j) e E 
\Sij\<8^ V(i, j), (j,i) e E 

tan(-0$)X(Wij) < tan(0$)K(Wij) V(i,j) e E 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 

(20f) 


Notice that constraint (19e) is a convex second-order cone constraint, which is widely supported by 
industrial strength convex optimization tools (e.g., Gurobi [7], CPlex [9], Mosek [23]). 

The complete SOC relaxation of (AC-E-PF) is presented in Model 2 (SOC-E-PF). Constraints (20a) 
capture KCL and constraints (20b)-(20c) capture line power flow in the IT-space. Constraints (20d) 
strengthen the relaxation with voltage product second-order cone constraint, and constraints (20e)-(20f) 
capture the line power flow and phase angle difference operational constraints. 

Lemma 3.1. (SOC-E-PF) is a relaxation of (AC-E-PF). 

Proof. A reduction from the SDP relaxation was first observed in [20]. □ 


3.2 The Network Flow Relaxation (NF) 

The inspiration for the network flow relaxation is to produce a relaxation that is similar to the linear flow 
models that are widely studied in operations research and computer science [1], This section proposes 
Model 3 (NF-E-PF) as such an analogue. The bulk of the operational constraints in this model are 
identical to (SOC-E-PF), however the convex non-linear thermal limit constraint (20e) is omitted in the 
interest of linearity. The key deference in this model is that the line flow constraints are simplified in to 
the linear expression (21b), and no longer require the Wij variables. The resulting model is essentially a 
traditional network flow [1] with two additional constraints to capture the line losses (21b) and the phase 
angle differences (21c)-(21d). Observe that, when there is no line charging, constraints (21b) simply 
express that the line losses must be nonnegative, i.e. power cannot be created on a line. 

In the rest of this section, we demonstrate that (NF-E-PF) is a linear relaxation of (SOC-E-PF). 
Through a series of deductions we will show that (21b) is simply a weaker version of (20b)-(20d), 
demonstrating the relaxation property between these two models. We begin by observing that the line 
losses in (SOC-E-PF) are given by, 


■ Sh = Y* 


Wi 

\Tn\ 2 


w y 

J-i* 


W* 


Wi - i-z- 


Wi 

IT, | 2 


Wi 


Lemma 3.2. (19e) ensures Wi/\Tij\ 2 — Wij/T)* — W*/Tij + Wj ^ 0. 


(22a) 


6 








Model 3 The Network Flow Relaxation of the Extended AC Power Flow (NF-E-PF) 


variables: 

Sfe(Sf\Sf u ) Vi e N 
Wtedvl) 2 ,^?) 2 ) Vi e N 
Sye^.Sg) V(i,j),(j,i)eE 

subject to: 

S?-St-(YS)*Wi= £ S ij+ £ Sji Vi e N (21a) 

( i,j)eE {j,i)sE 

Sij + S jt > -Ef + W^j (i,j) e E (21b) 

tan(-fl^)» (z*T* ((V,* - - S«)) < 

3 (ztjT* (Jy* - i b f) ^2 - Sii)) V(i, j) e E (21c) 

W®«)» ((*b - - %)) > 

3 (z*-!)* ((V,* - ^ - %)) V(i,j) g E (21d) 


Proof. First, observe that (19e) is equivalent to, 


W u 


rp* 

ij 


Wi 


which is equivalent to the following, 


SR 


Wi 

\Ta\ 2 " 3 


W, - A 


w± 

rp* 

*3 


Wi 

rp* 

ij 


+ A( 
SR 


(Wi 


PL 

\ T" 1 * 

\ Jm ij 

Wfl 

i-3 


Wi 1xr 

" \TiA 2 j 

Wi 

‘ rnlF 11 ''- 3 


w» 

rp* 


H SR 


W, 

ij 


_ 13_ 

rp* 


< 


Wi 

\Ti i \ 2,rj 


Wi - A 


Wi 

rp* 

ij 


Noticing that Wi, Wj, A(-) , |7\,j 2 are all positive, we can deduce the following bounds, 


Wi 


WA 2 3 


Wi sS 3R 


Wi 

TS 

ij 


_A? 

rp* 


Wi 


Now observe the equivalence, 

Wi 

frA 2 


Wi 

r 

ij 


_ 

rjp* 


W* 


+ Wi 


Wi 

\TiA 2 


+ Wi - 2 SR 


Wi 

li 

ij 


_ % 3 _ 

rjp* 


(23a) 

(24a) 

(24b) 

(24c) 

(25a) 

(26a) 


We want to determine the smallest possible value of (26a). Given that Wi,Wj are strictly positive, the 
largest possible value of SR(W'V f - /T*A will minimize (26a). Utilizing (25a) as a bound for this value, we 
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Model 4 The Copper Plate Relaxation of the Extended AC Power Flow (CP-E-PF) 


variables: 


S°e(S? l ,S? u ) Vi e N 
Wi e (K*) 2 ,«) 2 ) VzGiV 

subject to: 

Y i sf-T i s i-T i (y')*w i > 2 - {b \ 

ieN ieN ieN (i,j)sE 



(30a) 


have that smallest possible value of (26a) no smaller than, 


which factors to, 



Given that the square of any expression is a positive number, the result follows. 
Theorem 3.3. (21b) is a valid redundant constraint in (SOC-E-PF) when Zij ^ 0. 
Proof. Zij ^ 0 •«=> Y* ^ 0, combining this fact with Lemma 3.2 yields, 


(27a) 


(28a) 

□ 


Wi 

I ?o ] 2 


Wo 

rp * 

ij 


W* 

ij 


+ Wj)> 0 


The result follows from combining (29a) with the definition of line loss (22a). 
Lemma 3.4. (21c)-(21d) are valid redundant constraints in (SOC-E-PF). 


(29a) 

□ 


Proof. Begin by factoring (20b) similarly to (10). Then it is easy to observe that (21c)-(21d) can be 
derived by combining (20f)-(20b) in (SOC-E-PF). □ 


Corollary 3.5. (NF-E-PF) is a relaxation of (SOC-E-PF). 


3.3 The Copper Plate Relaxation (CP) 

The inspiration for this relaxation is to develop a version of the classic copper plate approximations used 
in power system analysis, which simply state that supply and demand should be balanced throughout 
the network. This section proposes Model 4 (CP-E-PF) as the relaxation version of a copper plate 
approximation. All of the operational constraints in this model are identical to (NF-E-PF), however the 
constraints and variables relating line flows (Sij) have been eliminated. Additionally, all of the KCL 
constraints have been combined into one simple linear expression (30a). A key benefit of this relaxation 
is that it is incredibly simple and scalable, while still capturing some properties of the lines. In the rest 
of this section we demonstrate that (CP-E-PF) is a relaxation of (NF-E-PF). 

Theorem 3.6. (30a) is a valid redundant constraint in (NF-E-PF). 











Proof. Due to the conjunctive nature of mathematical programs, combining constraints yields redundant 
constraints. As a first step, all of the KCL constraints (20a) are combined together. 


£ Sf - Sf - {Yf)*Wi 

ieN \ 


£ £ Sn 

(i,j)eE ( 


£s?-£s?-£(i?)*wi = £ 

ieN ieN ieN ieN 


£ Sy 

pi,i)eE 


£ Si* 

U,i)sE ) 


£ sf - £ sf - £ (Yi s )* w i = £ Sv + Sji 

ieN ieN ieN ( i,j)eE 


(31a) 

(31b) 

(31c) 


Observe that line loss constraints (21b) can be used to eliminate the S variables from this expression 
entirely, yielding, 


£sf-£s?-£(ir)*w;> £ -i 

ieN ieN ieN ( i,j)eE 


Wi 
\Tij \ 2 


w. 


(32a) 


which demonstrates the result. 


□ 


Corollary 3.7. (CP-E-PF) is a relaxation of (NF-E-PF). 

Note (CP-E-PF) includes all of the assumptions of (NF-E-PF) and hence can only be applied when 


4 Comparison of Linear Relaxations 

One of the key benefits of both (NF-E-PF) and (CP-E-PF) is that they are linear relaxations of (AC- 
E-PF). Although linear approximations of (AC-E-PF) are quite common [21, 5, 18], linear relaxations of 
(AC-E-PF) are, to the best of our knowledge, limited to two previous works [22, 4], which we review in 
detail. 

There are two key differences of this work and [4]. First, [4] focuses on only the rectangular real 
number formulation of (AC-E-PF) and all of the proofs follow from that formulation. This work uses 
complex numbers directly for the proofs, hence the results apply to any real number realization of these 
models. Second, the primary goal of [4] is to propose valid linear inequalities of non-convex constraints 
for iteratively strengthening relaxations of the (AC-E-PF), while this work focuses on building a intuitive 
static model. Indeed, the results of [4] may be used in consort with this work, to iteratively strengthen 
the linear relaxations proposed herein. 

The most closely related work to this is [22], which proposes a relaxation along the lines of (NF-E-PF). 
The key difference is that instead of constraints (21b), [22] uses the following valid equalities, 

• %))- ■ IF;) V(i,j) G £ (33a) 

»(*y (Sij - Sji)) = (Wf + - Wj} 6 E (33b) 

Because these constraints are equalities (not inequalities), it is difficult to develop a strong theoretical 
connection to the models proposed herein. However, as will be demonstrated in an example, the linear 
relaxation proposed in [22] can generate power on the lines, which is a significant disadvantage over 
the models proposed here. In the rest of of the paper TH will be used to denote the linear relaxation 
proposed in [22]. 

4.1 An Illustrative Example 

This section illustrates all of the power flow relaxations considered here on the 3-bus network from [14], 
which has proven to be an excellent test case for power flow relaxations. This system is depicted in Figure 
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Figure 1: 3-Bus Example Network Diagram. 
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Table 1: Three-Bus System Network Data (100 MVA Base). 


1 and the associated network parameters are given in Table 1. This network is designed to have very 
few binding constraints. Hence, the generator and line limits are set to large non-binding values, except 
for the thermal limit constraint on the line between buses 2 and 3, which is set to 50 MVA. In addition 
to its base configuration, we also consider this network with reduced phase angle difference bounds of 
18 degrees. The nonlinear global optimization solver Couenne [3] is used to find the globally optimal 
solution to (AC-E-PF) and the optimally gap between this solution and a relaxation is computed using 
the formula, 

Heuristic — Relaxation 
Heuristic 

Hence, the smaller the optimality gap the better the relaxation. 

Table 2 summarizes the results. As expected, in both cases the SOC relaxation has the smallest 
optimality gap. Further relaxing the model to NF or CP, increases the gap by about 1.5%. In this 
example the NF and CP produce the same optimally gap, however that is not the case in general. The 
TH relaxation has an optimality gap of over 80% indicating that it is significantly weaker than the other 
linear relaxations considered here. 


5 Conclusion 

This paper has developed two intuitive linear relaxations of AC power flow for transmission systems, one 
based on network flows (NF) and another inspired by copper plate approximations (CP). Although it 
was shown that both models are weaker than the established nonlinear SOC relaxation, their linearity 
brings significant performance and scalability benefits. Combining these linear relaxations with the 
established nonlinear relaxations (SDP [2], QC [8], SOC [10]) provides a rich variety of tradeoffs between 
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Table 2: AC-OPF Bounds using Relaxations on the 3-Bus Case. 


the relaxation accuracy and scalability. Given that all of these power flow relaxations span a wide range 
of mathematical programs and associated solution methods, the natural frontier for future work is to 
conduct a detailed numerical study to better understand the time and quality tradeoff of each relaxation. 
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A Real Number Formulations 

In the interest of clean proofs, this document focuses on a complex number representation of the relax¬ 
ations. In practice however, the a real number representation is often needed for implementation. This 
section presents the real number version of (NF-E-PF) and (CP-E-PF) to aid the implementations of 
these models. 

First we develop a real number short-hand for ZijTij as follows, 

K ( Z^Tij ) = tz^j = rijtj(j - Xijtfj ( i,j ) e E (34a) 

A ( Z^Tij) = tzjj = njtij + Xijt* (■ i,j ) e E (34b) 

Notice that the superscripts R, I are used for the real and imaginary parts of a complex number respec¬ 
tively. Expanding (10) using this short-hand into real and imaginary components K(V)V)*),Sy(V)V)*) we 

have, 

&(ViV*) = - - tzlj |r “ tz ijPii ~ (b J) e E (35a) 

A(UU/) = - tzf- ~ tz ij9ij + tzfjPij (■ i,j ) e E (35b) 

both of which are used to implement the phase angle difference constraints in any model with v,p,q 
variables. The complete real number implementation of (NF-E-PF) and (CP-E-PF) are presented in 
Model 5 and Model 6 respectively. 
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Model 5 The Network Flow Relaxation of the Extended AC Power Flow (NF-E-PF) in Real Numbers. 


variables: 

pf e (pf ,pf“) Vie N 
q?e(q? l ,q? u ) Vi e N 
w i e (K Z ) 2 >«) 2 ) VieN 
Pij e Sij) V(i,j), (, j, i) e E 
( «".*“) V(i,j),(j,i)eE 
subject to: 

Pi - Pi - 9i w i = XI Pi i + X Pji V * e N 

( i,j)eE (j,i)eE 

q 9 i - q? + b s z Wi = X Vij + X Va Vie N 

(i,j)eE (j,i)tE 

Pij T Pji ^ 0 (t, j ) £ E 

( lij + Qji 5 s (bi) e E 

tan(-0 A ) ~ **«P« - A 

( f, ~ tz ij (y')) £T - tz ij9ij + tz ijPij ('•./) 6 ^ 

tan(6» A ) - tzfj |r- - tzffpij - tz^q^j ^ 

(■-y)) - tz ij9ij + tz ijPij (b J) e ^ 


(36a) 

(36b) 

(36c) 

(36d) 


(36e) 


(36f) 


Model 6 The Copper Plate Relaxation of the Extended AC Power Flow (CP-E-PF) in Real Numbers. 


variables: 

pf e (pf ,pf“) Vi e N 
q 9 £ (qf,qf u ) VieN 
Wi e((vlf,(vT) 2 ) VieN 

subject to: 

Xpf - X ^ _ X u,i ^ 0 ( 37a ) 

ieiV zeiV ieJV 

X q i - X + X > X (tT + (37b) 

ieiV i€N ieN ( i,j)eE V b / 
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